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Abstract 

A vortex line, shaped by a zigzag of pinning centers, is described here through a three-dimensional 
unit cell containing two pinning centers positioned symmetrically with respect to its center. The 
unit cell is a cube of side L = 12£, the pinning centers are insulating spheres of radius R, taken 
within the range 0.2£ to 3.0£, £ being the coherence length. We calculate the free energy density 
of these systems in the framework of the Ginzburg-Landau theory. 
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I. INTRODUCTION 



Superconductors have many kinds of imperfections, that is, internal regions where Cooper 
pairs either don't exist at all or exist purely as a fluctuation effect. Thus the macroscopic 
wave function, describing the collective state, vanishes abruptly or asymptotically inside such 
regions. In presence of an external magnetic field vortices arise inside the superconductor 
and are strongly attracted to such imperfections, also called pinning centers. For this reason 
pinning centers have been extensively studied in the past in many ways, including artificially 
made ones, such as columnar defects 1 -, antidotal and micro holes 4 . They are interesting 
because they bring clear-cut questions about the interaction between vortices and pinning 
centers^, such as how local misalignment really occurs inside the superconductor. Vortices 
lines should be aligned to the applied field but the presence of strong attraction to a pinning 
center can change this locally. As a result of competing energetic demands new interesting 
phenomena can take place in vortex Physics, such as the one considered here. A vortex line 
in the absence of pinning centers is aligned along the magnetic induction direction, hereafter 
called z-axis. The presence of a zigzag of pinning centers forces the vortex line to bend and 
acquire this shape, resulting into local misalignment, though it remains oriented along the 
magnetic induction. 

Pinning forces act on the vortex core whose radius is given by the coherence length £. 
The interaction of pinning centers with vortices has been studied using several approaches 7 . 
From the point of view of the Ginzburg-Landau theory, pinning may be caused by spatial 
fluctuations of the critical temperature^, T c (x), or of the mean free-path^ that changes the 
coefficient in front of the gradient term, £(x) 2 |(V — ^A)A\ 2 . The interaction between a 
vortex line and a pinning center has been considered by many authors in the context of 
the Ginzburg-Landau theory*^. The number of vortex that can be trapped by a defect 
is an interesting problem. In case of a columnar defect the saturation number has been 
determined long ago by Mkrtchyan and Shmidlji 1 -: n s = A/2£, though this formula has to 
change near the upper critical fielcU 2 -. Recently the saturation number has been discussed 
for three dimensional cavities^. 

The system under investigation here consists of a vortex line trapped by a zigzag of very 
large pinning centers, namely, spherical insulating cavities, with radius taken to vary from 
0.2£ to 3.0£. For these pinning centers the boundary-value problem has to be taken into 
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account, as usually treated since the de Gennes boundary conditioniamust be satisfied at the 
cavity surface. We chose to describe this system through a unit cell with two cavities inside, 
that rotate freely around its center, producing for each angle a distinct zigzag arrangement. 
In this paper we analyze numerically the angular dependence of the Helmholtz free energy 
T c {9, R) and propose an expression for it. 

This paper is organized as follows. In section |H] we present the model for a 3-D super- 
conducting media with the cavities. In section IIII1 we discuss out theoretical approach. In 
sections IIVI and IS we show the results obtained through numerical simulations. In section 
IVII we summarize the main results of the work. 

II. THE MODEL 

The system studied here consists of a cubic unit cell with size L equal to 12£, and 
two insulating cavities of radius R, as specified by the figure d The line segment joining 
the cavities makes an angle 9 ranging from 0° to 180°, taken here in increments of 3°. 
Notice that because of symmetry 9 varying from 0° to 90° is sufficient to obtain all possible 
configurations. 
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FIG. 1: The plane of rotation of the cavities inside the unit cell. The distance D between the two 
cavities is equal to L/2. The superconductor fills the remain of the unit cell, shown here as a gray 
region. 

Figure^ shows the plane of rotation (plane y) and parameters associated to the unit cell. 
While 9 and R are freely changed, D, the distance between the center of the cavities, is 
fixed and equal to L/2. Thus the present model features a distance between the two cavities 
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independent of the cavity radius. 

Our numerical treatment demands a mesh grid to describe the unit cell. The length of 
the cube side is 12£, and we choose the number of mesh points along a given direction, 
P, to be 19. This choice implies that distance between two consecutive mesh points, a, be 
equal to 2£/3. The value of a must be smaller than the coherence length £, which is the 
minimum physical scale of the Ginzburg-Landau theory. The number of grid points inside 
the cavity should be large enough to describe it. This number is obtained from the ratio 
between the two volumes, of the cavity and of the cube. Since the entire cube has P 3 points, 
the number of grid points inside the cavity is |7r(-^) 3 (l + -|) 3 . For instance, a cavity of 
radius R = 1.0£ has approximately 16 points, whereas the R = 2.0£ sphere has 8 times more 
points. Obviously in the limit of a cavity with radius smaller thant the mesh distance, that 
is R < a = 2/3£, there will be just one point in the mesh describing the cavity. In this case 
the order parameter does not vanish inside the defect but just undergoes a drop on its value. 

III. THEORETICAL APPROACH 

We start our considerations for the energy density functional of the Ginzburg-Landau 
theory^, expressing it in units of the critical field energy density^ 3 -, H^/Air. The periodicity 
of the problem requires a search of the free energy minimum for a fixed integer, the number 
of vortices inside the unit cell. 

This integer also fixes the magnetic induction B{x), which is the average of the local field 
taken over the unit cell volume, 



The magnetic induction is completely determined by the vorticity of the system because 
the overall current circulation vanishes inside the unit cell. The relationship between the 
vorticity of the system u(f) and the magnetic induction B(x) in reduced units is, 



where V = n x x + n y y + n z z is the vorticity in an arbitrary direction. In the present paper 
we consider v = z and present some results concerning the V = 2z case. The parameter 
k = A/£ is the dimensionless Ginzburg-Landau parameter and A is the penetration depth. 




(1) 




(2) 



4 



Here we consider the no magnetic shielding limit. The field penetrates in the superconductor 
with no Meissner-Ochsenfeld effect. In this regime h(x) = V x A(x) = B. This situation 
can be viewed as a large k limit. In reduced units the free energy density is normalized by 
the critical field density, H^/4ir, and the order parameter density is dimensionless, varying 
between and 1. 

2ni 2 



-IAI 2 



+ §|A| 4 , (3) 



The function t(x) is a step-like function used to specify the cavities in this approach^ 3 -. 
Explicitly we have t(x) = t\{x)t2{x) and 

^ = 1 " i + e (l*-*l/*)" ' (4) 

where Tj is equal to inside and 1 outside the ith cavity. The above explicit representation 
of the t function is necessary for computational reasons and for accuracy we take that 
N = 8. In the limit N — > oo, the function r tends to the well-known Heaviside function, 



R J \ R 

Since we are in the no shielding limit, the vector potential A(x) is determined from Eqs. 
[I] 121 and the condition of magnetic flux quantization inside the unit cell. The vector potential 
does not participate in the minimization process of the free energy density that only takes 
into account the real and imaginary parts of the order parameter. The free energy density 
contains two terms. The first term is the condensation energy density, —t(x) |A| 2 + | |A| , 
which in case of no vortices (|A| 2 = 1) and no cavity (r = 1) has the value -0.5. The 
presence of a cavity raises the energy since inside it the density vanishes (|A| 2 = 0). And 



the second term, the kinetic energy density, r(x)£ 2 (V — ^A)A 
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Notice that there is 

kinetic energy in case of no vortices but with a cavity. At the insulating-superconducting 
interface r changes from 1 to and this causes a bending of the order parameter, which has 
some kinetic energy cost. 

The most significant advantage of the present method, is that the free energy functional, 
Eq. El contains the appropriate boundaries conditions to the problem. This removes the 
necessity of solving the theory in two independent regions and later applying the Neumman 
boundary conditions. Besides the present method easily applies to internal regions of any 
shape, not just spherical, and finds its solution for the given normal-superconductor interface. 
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IV. RESULTS 



For a given pinning arrangement, which means fixed values of 9 and R, and the condition 
of one vortex in the unit cell, we carry the minimization of the free energy. Here we present 
the results of several independent simulations obtained for the pairs (9,R), ranging from 
(0°, 0.2£) to (180°, 3.0£), in increments of 3° for 9, and of 0.2£ for R. For each simulation we 
initialize the order parameter in a random way and the minimization procedure is carried for 
each temperature^. Then the temperature is lowered until a convergence criteria is reached, 
which means that the free energy has become stable. For each simulation we implement at 
least 1200 Monte Carlo visits per mesh point in a Metropolis algorithm. 



Our main results are shown in the two main curves of figure EJ The figure 2(a) shows 
several curves of the free energy density versus the angle 9, each curve associated to a different 
R value. Similarly the figure |2(b)| shows several free energy density versus R curves, each 
corresponding to a distinct value of 9. The figure 12(a)] shows that mirror symmetry 9 tt — 9 
holds as expected, since the two cavities inside the unit cell are equivalent. The free energy 
T c {9, R) is an even function in 9. 

In our previous wor k 16 i 17 we have found the remarkable property that cavities inside a 
superconductor can lower its energy as compared to the cavity-free superconductor. This 
effect can be verified here in figure |2(b)[ which display a set of points lying under the 
energy threshold of —0.434, the free energy !F C (9,R = 0), of the system without cavities, 
approximated by T c {9 — 0, R = 0.2), as previously discussed. In summary we found here 
several pinning configurations, each described by the pair of values (9,R), that have lower 
energy than the cavity free system. Table exhibits these pairs (9, R) of lower energy. 



TABLE I: Pair of values (0, R) that establishing a configuration with lower energy than the cavity- 
free superconductor. 
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The curve of figure |2(b)| shows monotonic growth^ for the free energy, but only for R 
equal or larger than 1.2£. It also shows a local maximum for the free energy at R = 1.0£. 
In fact the behavior is distinct for the two regimes, R > 1.0£ and R < 1.0£. To understand 
the R > 1.0£ case lut us hold the angle fixed, and vary the cavity radius. One finds that the 
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FIG. 2: Dependence of the free energy, J- c , with the pinning center radius R and the angle 9. 



2(a) Variation of T c as a function of the angle 9 in the range 0° to 180°, data points obtained for 
increments of 3°. The radius R varies from 0.2£ to 3.0£ and for each increment of 0.2£ results in a 
distinct curve, all plotted in ascendant order from bottom to top. 



2(b) Variation of J- c as a function of the pinning center radius for a specific value of 9. The radius 
varies from 0.2£ to 3.0£ with an increment of 0.2£. Distinct curve correspond to different 9, equal 
to 0°, 9°, 18°, 27°, 36°, 45°, 54°, 63° and 72° in the ascendant order form bottom to top. 

free energy density is given by a constand density times the volume of the unit cell, Vq = L 3 , 
removed of the non-superconductor volume of the cavities, V c = 8/3nR 3 . This result is 
only approximately valid since the curvature of the order parameter near the pinning sphere 
surface causes an increase in the kinetic energy, an effect that becomes more pronounced for 
large spheres. Thus we have that for R > 1.0£, 

/ 8irR 3 \ 



fixed, R) OC T c f 1 



3L 3 



(5) 



where T c f = ^(0°, 0.2£) = —0.434 is the cavity-free energy. The monotonic growth seen in 
figure |2(b)| is well described by the cubic R dependence of Eq. El 

The dependence with the angle 6 exhibited in figure |2(a)| may be described by the ex- 
pression 



!F C (9, Rfixed) OC fk(Rfixed) SUl 2 6 



(6) 
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where the function fk(R) incorporate all the kinetics effects produced by the presence of 
the cavities with radius R. In figure |2(aj| fk(R) provides the amplitude of oscillation for a 
fixed R. As we have mentioned before, the kinetic effects become more pronounceable as 
the cavity radius increase. We assume that fk is a linear function of R. 

We add Equations and El together to obtain a general free energy expression to describe 
arbitrary angles of rotation and radii bigger than the coherence length: 

T c {9, R) = T c} f k (R) sin 2 9. (7) 



V , 

For all radii, the configurations of minimum and maximum energies are obtained for 0° 
and 90°, respectively. The maximum at 9 = 90° has a smaller superconducting volume as 
compared to the minimum at 9 = 0° configuration. This is easy to understand because 
the 0° configuration has the two spheres aligned along the z-axis and both overlap with the 
vortex line, whereas the 90° configuration only one overlaps the vortex line. The other one is 
free in space thus taking away space that could be otherwise superconducting. This makes 
the 90° configuration closer to the normal state than the 0° one. At some intermediate angle 
between 0° and 90° a depinning transition takes place, although it is not noticeable in both 
figures |2(a)| and |2(b)| This transition has been studied in ReiM. The function fk(R) is 
easily obtained by taking its difference at extreme angles, T c {9 max , R) — T c {9 min , R), where 
9 max = 90° and 9 min = 0°. This difference is shown in figure El as a function of R. Thus 
Eq. the major conclusion of this paper, gives a good description of the free energy for all 
(6, R) pairs. 

In the limit that R — > 0, the free energy should converge to the cavity free superconductor, 
J-'cf, whose energy only depends on the unit cell vorticity. Thus the term fk should vanish 
in this limit R — > so that the angular dependence sin 2 9 disappears from the free energy: 

Kmfk(R) -0. (8) 

We determine the fk{R) term fitting the curve of the figure El with the best linear function. 
In this way, the function found is: 

f k (R) = -0.00358 + 0.00962.R. (9) 

The linear dependence of the Eq. 01 with the radius R reflects the importance of the kinetic 
energy, £ 2 (V — ^A)A , which describes the bending of the order parameter at the surface 
of the cavities, that becomes more important for large cavities. 
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FIG. 3: Amplitude of oscillation. JF c (90°,i2) - f c (0° , R) 

Substituting all the terms in Eq. [7|by the terms obtained in the fitting and the parameters 
of the system, the energy dependence is expressed by following equation: 

T c {9 y R) = F cf (l - + (0.00962/? - 0.00358) sin 2 9. (10) 

The negative constant —0.00358 is in conflict with the condition expressed by Eq. [HI In fact 
in the region R < 1.0£, Eq. El does not apply because of mesh effects. Besides for such small 
R the present Ginzburg-Landau approach is not applicable, as previously discussed. 



V. MANY VORTICES NEAR THE ZIGZAG OF CAVITIES 

So far we have described the case of just one vortex near a zigzag of cavities, v = lz. In 
this section we briefly comment on the general case of many vortices along the z direction, 
V = n z z. As n z increases the upper critical field is approached and beyond the upper critical 
field there is also a surface superconductivity state at the surface of the cavities. The radius 
of the cavities R and the angle of rotation 9 sets new geometric configurations in the unit 
cell that lead to multiple trapping and giant vortex states in the superconductor—. A simple 
example of such configurations is shown in figure 0J which displays two vortices in the unit 
cell (n z = 2 case) for R = 1.8£ at the two extreme angles, 9 = 0°, and 9 = 90°. For 
9 = 0° the system presents strong competition between the vortex-vortex repulsion and the 
vortex-cavity attraction. While one of the vortices is barely trapped by the cavities the 
other vortex is not pinned. The rotation of the cavities to the 9 = 90° configuration yields 
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(a) (b) 
FIG. 4: Iso-surfaces of the order parameter |A(x)| 2 inside the unit cell, a cube with side L = 12£ 



containing cavities (in red) with radius R = 1.8£. The figure 4(a) corresponds to 9 = 0° and shows 
just one vortex line trapped by the defect, whereas the other one is not pinned and repelled due 
to vortex- vortex repulsion. In the figure |4~(b)| the cavities form an angle = 90° and each vortex 
is trapped by an independent cavity. 

a state of two vortices each one trapped by a distinct cavity. For values of n z bigger than 2 
more complex situations are possible and will be studied elsewhere. 

VI. CONCLUSIONS 

We have carried here a Ginzburg-Landau theory study of a vortex line near a zigzag of 
pinning centers. The pinning centers are insulating spherical cavities and their arrangement 
is well described by a unit cell containing two of them. Inside the unit cell the rotation 
of the two cavities around the center produces a continuous of zigzag arrangements whose 
effects on the vortex line are here discussed. We observe distinct behavior for cavity radius 
above and below 1.0£. Sweeping the angle of rotation of the cavities around its center, 9, 
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makes the zigzag more pronounced to the point that the vortex line decouples from it above 
the critical angle 9 C , as previously found 13 . Below this angle the vortex line is pinned by 
both cavities and above by just one. Here we have determined that the free energy density 
has a simple analytical dependence, expressed by the Eq. [7| 
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